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Investigated herein is the buckling of nanobeams based on a nonlocal Timoshenko beam model 
by the method of initial values within the framework of nonlocal elasticity. Since the nonlocal 
Timoshenko beam theory is of higher order than the nonlocal Euler—Bernoulli beam theory, it is 
known to be superior in predicting the small-scale effect. The buckling determinants and critical 
loads for bars with various kinds of supports are presented. The Carry-Over matrix (Transverse 
Matrix) is presented and the priorities of the method of initial values are depicted. To the best of 
the researchers’ knowledge, this is the first work that investigates the buckling of nonlocal 
Timoshenko beam with the method of initial values. 
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1. Introduction 


The Timoshenko beam theory allows the rotation between the cross-section and the 
bending line. This deformation is related to shear rotation, which is not considered in 
a Bernoulli beam. As a result, shear strains and stresses are removed from the 
Bernoulli beam theory. Therefore, Bernoulli beam is powerful. The Timoshenko 
beam theory was presented in the early part of the 20th century. Nowadays, the 
buckling of Timoshenko beams is investigated by experimental and numerical 
methods.' Wu et al.” investigate the free vibration and elastic buckling of sandwich 
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beams with a stiff core and functionally graded carbon nanotube reinforced 
composite (FG-CNTRC) face sheets within the framework of Timoshenko beam 
theory. 

Today, interest in nanotechnology is increasing strongly, and many nanos- 
tructures are produced. A nanostructure is a small object of intermediate size be- 
tween molecular and microscopic (micrometer-sized) structures. Structural elements 
such as beams, plates, and membranes in micro- or nanolength scale are widely used 
as components in micro- or nanoelectromechanical system (MEMS or NEMS) 
devices.*° To better understand the mechanical behavior of these structures; small- 
scale effects should be considered. The classical continuum theories are not capable of 
capturing the small-scale effects when dealing with nanoscale structures. In 1972, 
Eringen’® presented the nonlocal continuum mechanics to take into account the 
small-scale effect. According to this theory, the stress at a reference point in a body 
depends not only on the strains on that point, but also on strains on all other points 
of the body. Peddieson et al.° first applied the nonlocal continuum theory to the 
nanotechnology in which the static deformations of beam structures were obtained 
by using a simplified nonlocal beam model based on the nonlocal elasticity theory of 
Eringen.'° Later, many researchers examined the bending, buckling, vibration and 
wave propagation problems using the nonlocal elasticity theory.‘‘ 1° Wang and Hu”? 
studied flexural wave propagation in single-walled carbon nanotubes (SWCNTs) and 
presented that only the nonlocal elastic Timoshenko beam is able to predict the 
decrease in phase velocity when the wave number is so large that the microstructure 
of carbon nanotubes has a significant influence on the flexural wave dispersion. Wang 
et al.?' used the principle of virtual work to obtain the governing equations and the 
boundary conditions for the buckling of axially loaded rods/tubes based on the 
nonlocal Timoshenko beam theory. In this paper, the method of initial values 
(transfer matrix method) is used to obtain the critical loads of rods with various end 
conditions in the frame of nonlocal elasticity. This method gives the values of the 
displacements and stress resultants throughout the rod once the initial displacements 
and initial stress resultants are known. A priority of the method is that the high 
degree of statical indeterminacy adds no additional hardship to the solution of the 


problem.*?? 


2. The Physical Meaning of the Nonlocality Parameter and 
a Short Review of Nonlocal Elastic Models in Modeling 
of Carbon Nanotubes and Graphenes 


Nonlocality parameter can have different physical meanings: a distance between 
particles in discrete structures, the dimension of the grain or a cell, a characteristic 
radius of correlation or action — at — a distance forces, etc. However, we shall 
always assume that the scale parameters are small in comparison with the dimen- 
sions of the body.”° According to the aforementioned studies on the basis of Eringen’s 
nonlocal elasticity theory, a consensus has been reached that the refined nonlocal 
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models are superior to their local counterparts. However, since the theory also 
belongs to the category of continuum mechanics, verification of nonlocal continuum 
models is essential to identify the magnitude of the small-scale parameter y which is 
the key revealing scale effect in models.** The nonlocality parameter revealing the 
small-scale effect on the responses of structures of nanosize. In recent years, a large 
number of studies have been conducted to determine the nonlocality parameter.””° 
Generally, a conservative estimate of the nonlocality parameter is y < 2nm for a 
SWCNT for wave propagation.” Wang'® indicated that nonlocal elasticity theory 
should be used for an accurate prediction of wave propagation in carbon nanotubes 
(CNTs), and estimated a range of the small-scale parameter. Khademolhosseini”® 
examined the vibration of CNTs with nonlocal elasticity formulation and molecular 
dynamics simulation, showing that both models are in great agreement for non- 
locality parameter (y) is 0.185nm. Liang and Han”? show that nonlocal scale pa- 
rameter is independent of the length of CNTs, and dependent on the radius of CNTs. 
In Han and Weng’s paper,”° the length-dependent stiffness is revealed by the non- 
local elasticity and verified through molecular simulation results. The value of the 
scale coefficient in the nonlocal theory is recommended to be about 0.7 nm for the 
application of the nonlocal theory in analysis of CNTs. Arash and Wang?’ present a 
review of the application of nonlocal elastic models in modeling of CNTs and gra- 
phenes. Explicit expressions for the critical buckling load and strain were derived for 
axially loaded rods/tubes based on the nonlocal Timoshenko beam theory.?!°°°? 
Wang et al.?!°°?! developed nonlocal elastic beam and shell models to investigate 
the small-scale effect on buckling analysis of CNTs under compression. 

The nonlocal Timoshenko beam model was implemented**~*° to investigate an 
axial buckling of SWCNTs embedded in an elastic medium, and the influence of the 
small-scale effect and stiffness of the surrounding medium was investigated. Zhang 
et al.°°" presented an assessment of nonlocal beam and shell models in the predic- 
tion of critical buckling strains of axially loaded SWCNTs by use of molecular 
dynamic (MD) simulation results for SWCNTs. 


3. Basic Equations of Nonlocal Timoshenko Beam and the Method 
of Initial Values 


Basic equations of classical Timoshenko beam are given by”! 


dv 1 
= T-Q 1 
dz «GF , (1) 
dl M 
dM 
dT d2v 
ar = aa (4) 
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where z is the longitudinal coordinate, v is the transverse displacement, 1) is 
the rotation due to bending, M is bending moment, T is the shear force, P is an 
axial compressive load, G is the shear modulus, F is the elasticity modulus, F' is 
the cross-sectional area of the beam and IJ is the second moment of area. The 
shear correction factor («) is given for circular and rectangular cross-sections as 
follows:°° 


6(1 + v) , _ 0+») 
T+6v ’ ~ 124+11p" 


(5) 


The following relation exists between local and nonlocal stresses. 
Ononlne( 2) = *y” @ nonitoe| 2) = Ging (2), (6) 


where + is nonlocality parameter. If the two sides of Eq. (6) are multiplied by ydF 
and integrated over the cross-section, the following relation is obtained: 


[ cwncl2uaF — 7? | Grane udF = f oe(2uaF (7) 
F F F 


By modifying this equation, one obtains 
Mosnioe(2) ~ 77M" sonloc(2) = Myoc(2). (8) 


By using Eq. (8), Eqs. (1)—(4) take the following form: 


dv 1 
dz ~ FGe . (9) 
dQ P— FG 
dz FGk(7?P—E]D+EIP ay 
dM 
aT FGKP (12) 


= M 
dz FGk(y?P—-EI)+EIP 


The system of equations given above is the basic equations of nonlocal Timoshenko 
beam subjected to an axial compressive load. 


(a) Equations (9)—(12) revert to the basic equations of nonlocal Euler beam in the 
limit « goes to infinity. 

(b) Equations (9)—(12) revert to the basic equations of classical Timoshenko beam 
in the limit y goes to zero. 

(c) Equations (9)—(12) revert to the basic equations of classical Euler beam in the 
limit y goes to zero and « goes to infinity. 
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The system of equations given by (9)-(12) can be written in matrix form as 


follows: 
0 -1 0 : 
FGk 
UV 
P— FGk 
4|9)_|° © perpen. ° 
a 0 0 0 1 
= FGP 
0 0 a 0 
EIP + FG(P7? — ED)« 
where 
Opto 0 
2 ; P— FG«k 
go) A. as EIP + FG(Py? — ED 
a 0 0 0 
FGPk 
0 0 . 


EIP + FG(P7? — EI)« 


The initial value problem can be written in the following form: 


v(0) v0 
Hag, wal (=| 2 
a M(0) Mo 
T(0) To 


Solution of initial value problem given by Eq. (15) is 
¥(z) = TM(z)o 


i. @ 
T 
1 
FG 
(14) 
1 
0 
(15) 
(16) 


(see Ref. 39), where TM(z) is called as carry-over matrix or transverse matrix. 


TM(z) is given below: 


1 cos(za) — 1 Elasin(za) — Pz(a?y? + 1) 
. P EIPo? 
0 (a2? + 1) sin(za) (a?2y? + 1)(cos(za) — 1) 
TM(z, a, El) = Ela Elo? 
0 O cos(za) a) 
a 
0 0 —asin(za) cos(za@) 
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where 


— VEVGViVP 
~ \/FGK(EI — 7?P) — EI P" 


(18) 


4. Examples 
4.1. Calculation of critical buckling load of rods/tubes with pinned ends 
In this case, two of the initial values are known (see Fig. 1). 
vj =0, My =0. (19) 
End conditions are 
o(L) =0, M(L)=0. (20) 


The remaining two initial values will be found with the help of end conditions. 
The solution of the problem can be written as follows with the aid of Eq. (16). 


v(2) = TM y2(2)Q% + TMya(z)To 


(21) 
O(z) = TMg9(z)Q% + TMo4(z)To (22) 
M(z) = TM39(z)Q + TM34(z)To (23) 
T(z) = TM42(2)% + TM ya (z)To (24) 


With the help of the end conditions, the buckling determinant is obtained as 


TM),(L) TMy(L) | Lsin(oL) _ 0sSaL=n7 (25) 
TMs)(L) TMz4(L) | Oo 7 


Critical buckling load of nonlocal Timoshenko beam (P) is obtained as follows (see 
Eq. (18)): 
w?EI 


P= = ; 
Dei t+ yr +1) 


(26) 


(a) The critical load given by Eq. (26) returns to the critical load of the classical 
Timoshenko beam (P,) in the limit 7 goes to zero. 


_ EIr? 
~~ 79 | Eln 
D+ ae 


/ El = 
L 
Se 


Fig. 1. Beam with pinned ends. 


P (27) 
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(b) The critical load given by Eq. (26) returns to the critical load of the classical 
Euler beam (P;) in the limit - goes to zero and « goes to infinity. 
mw? EI 
2= 73 
(c) The critical load given by Eq. (26) returns to the critical load of the nonlocal 
Euler beam (P3) in the limit « goes to infinity. 


(28) 


_ EL (29) 
3 yd 4 2° 
(see Ref. 21). 
4.2. Calculation of critical buckling load of rods/tubes with 
clamped-free ends 
In this case, two of the initial values are known (see Fig. 2). 
Uy =0, OQ =0. (30) 
End conditions are?! 
d 
To=-P—| = MOH. (31) 


dz z=L 


The remaining two initial values will be found with the help of end conditions. 
The solution of the problem can be written as follows with the aid of Eq. (16). 


v(z) = TM3(z)Mo + TMya(2)Tp, (32) 
Q(z) = TM3(z)Mo + TMo4(2)To, (33) 
M(z) = TM33(z)Mo + TM3,(2)To, (34) 
T(z) = TMy3(z)Mo + TMy4(2)Tp- (35) 
With the help of the end conditions, the buckling determinant is obtained as 
TMi3(L) — PTM'3(L)  TMyy(L) — PTM4(L) 
TM33(L) TM34(L) 
P(a*y? +1) cos(aL) _ _t 
5 ET Ss ol (36) 


Fig. 2. Beam with clamped-free ends. 
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Critical buckling load of nonlocal Timoshenko beam (P) is obtained as follows 
(see Eq. (18)) 
WEI 


P= — 
1? geqaz + “pr +4) 


(37) 


(a) The critical load given by Eq. (37) returns to the critical load of the classical 
Timoshenko beam (P,) in the limit 7 goes to zero. 


TRI 


f= Far : 
rH + 4p? 


(38) 


(b) The critical load given by Eq. (27) returns to the critical load of the classical 
Euler beam (P;) in the limit - goes to zero and « goes to infinity. 


_ WEI 


2 "4E2 - (39) 


(c) The critical load given by Eq. (27) returns to the critical load of the nonlocal 
Euler beam (P3) in the limit « goes to infinity. 
w?EI 


P, = ——.~——__-.. 40 
3 ry? + 4D? 0) 


4.3. Calculation of critical buckling load of rods/tubes with variable 
cross-section 


In the beam shown in Fig. 3, two of the initial values are known. 


Vo => 0, Mo = 0 (41) 
End conditions are 
vo(Ly + Ly) = 0, M (Ly + Ly) = 0. (42) 
ELF, EI,, F 
P 


Fig. 3. Variable cross-section beam. 
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The vector 7 at z= L, and z = Ly can be written as follows: 


Y(L1) = TM (Ly, oy, El, yp 
Y(L, + Ly) = TM (Lo, a2, Ely) y(L1) 


= TM(Lp, a2, Ely) TM (Ly, 01, Ely) yo. (43) 
The matrix TT'M is defined as 
ITM = TM (Lp, a2, Ely) TM (Lj, 0, El), (44) 
where 
1 cos(za;) — 1 EI, a; sin(za;) — Pz(y2a7 +1) 
P PEI, a? 
op (Pai +1) sin(zes) (y2aj + 1)(cos(zai) — 1) 
TM(z, a;, EI,) = EI,q; El,a? 
0 0 cos(za;) sizes) 
Oy; 
0 O a;(— sin(za;)) cos(za;) 
(45) 
The buckling determinant is obtained as 
(Ly + Ly) (ay sin(Qz D2) cos(Q, I) 
TIM TIM, _ + (aD) sin(a, I) cos(@ Dy)) ash (46) 
TTM3, TTM3, | — 1 ~ 
where the parameter a is as follows: 
(a) For nonlocal Timoshenko beam is 
____ VGvevP VF, i 
*  /GKE(EL, — 7?P) — PEI, 
(b) For nonlocal Euler beam is 
VGVPVF, 
Q; = ; (48) 
V-GR(?P — El) 
(c) For local Timoshenko beam is 
_ VGVaV/P YF; (49) 
*  ./-EI,(P — GeF) 
(d) For local Euler beam is 
P 
a; = vP : (50) 
EI, 
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The dimensionless critical loads are given in Figs. 4 and 5 for various values of the 
parameters. (a 2, W; H, B es 0 a, "% By €= 4, }; is the width of the 
cross-section and h; is the height of the cross-section.) 

(Porg: dimensionless critical load of local Euler beam, Ppy7: dimensionless critical 
load of local Timoshenko beam, Ppyg: dimensionless critical load of nonlocal Euler 


beam, Ppyr: dimensionless critical load of nonlocal Timoshenko beam). 


3 
ances, Po 
Ghs 
pero 
Por 
0.42 


0.40 


0.38 


0.36 


0.34 


0.00 0.02 0.04 0.06 0.08 0.10 $ 
Fig. 4. Dimensionless critical loads of variable cross-section beam. 
L3 
amet Pe 
Gr * 
0.145 
Poe 
0. 
0. Ppit 
0. 
0. 
0. 
0. 
i Ppne 
- Pont 


0.105 é 
0.02 0.04 0.06 0.08 0.1 


Fig. 5. Dimensionless critical loads of variable cross-section beam. 
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4.3.1. Special case 


The buckling determinant given in Eq. (46) reverts to the determinant in Eq. (25) 
for Ty = Dy» => L/2, F, = Fy and El, = Ely. 


4.4. Calculation of critical buckling load of Gerber beam 


Gerber beam is made by joining simple beams with hinges and removing the 
remaining freedoms by supports (see Fig. 6). 
In the Gerber beam shown in Fig. 6, the following values are known: 


v(0)=0, M(0)=0, M(Ly)=0, v(L) =0, A(L) =0. (51) 
That is 
0 v(L) v(Ly) 
0) =} | gaya] OO |, aya] OOF] Ge) 
(0) (1) T(L,) 


Using Eq. (16), the following relations can be written: 


j(L1~) = TM(L)).¥(0), (53a) 
i(L) = TM(L2).g(L1*) = TM(Ly).(y(L1-) + %), (53b) 
where 
0 
25 
=| 4 (54) 
0 


Using Eqs. (53a) and (53b), the following relations can be written: 
W(L) = TM(Ly).y(Ly~) + TM(Ly).t¢ 
= TM(L).9(0) + TM(Ly).a, (55) 
M(L,)=0, w(L)=0,) QA(L)=0. (56) 


Fig. 6. Gerber beam. 
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The following relations are obtained by the conditions given in Eq. (56). 
Ty(aEI sin(aL) — LP(a?y? + 1)) 


= } LQ = 0, 
o(L) Lp 72 ETP (0) =0 
Ty(a2y2 + 1)(cos(aL) — 1) 
+ A+ 2(0) = 0, 
aL) ober (0) 
Tp sin(L 
M(L;) = Tosin(Lia) _ 9 
a 
The buckling determinant can be written as follows: 
L L Elasin(La) — LP(a?y? + 1) 
. EIPa? 
Heel 4 , (a9? + 1)(cos(La) —1) | _ (L— Ly) sin(aLy) _ 0. 
Ela? a 
0 0 sin(L,q) 
a 


Critical buckling load of nonlocal Timoshenko beam is obtained as 


nw2EI 


Pyp = 


9(_m2EI 4, 727? ; 
Ly2( age + e+ 1) 


0.02 0.04 0.06 0.08 0.10 


Fig. 7. Dimensionless critical loads of Gerber beam. 


1950036-12 


(57) 


(58) 


(59) 


Buckling of Nonlocal Timoshenko Beam 


Critical buckling load of nonlocal Euler beam is 
nr EI 


= jj PB, = —————_.. 0 
Pyr tim NT my? +L? (6 ) 
Critical buckling load of local Timoshenko beam is 
nw EI 
Pur = (61) 


0.82 


0.80 


0.78 


0.02 0.04 0.06 0.08 0.10 


Fig. 8. Dimensionless critical loads of Gerber beam. 
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Fig. 9. Dimensionless critical loads of Gerber beam. 
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Critical buckling load of local Euler beam is 


neEI 


=r (62) 


Pur = lim Pxyp = 
7-0 


The dimensionless critical loads are given in Figs. 7-12 for various values of the 
parameters. (Pprp: critical buckling load of dimensionless local Euler beam, Pp: 
critical buckling load of dimensionless local Timoshenko beam, Ppyp: critical 


Ly 


Apa cr 
Gh Pope 


0.215 


0.210 


0.205 


0.200 


0.195 


0.190 


0.185 


0.180 


PpNE 


Ppt 


& 


p204 Ppnt 


0.02 0.04 0.06 0.08 0.10 


Fig. 10. Dimensionless critical loads of Gerber beam. 


PDLE 


PpNE 
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0.02 0.04 0.06 0.08 0.10 


Fig. 11. Dimensionless critical loads of Gerber beam. 
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can 


0.02 0.04 0.06 0.08 0.107 


Fig. 12. Dimensionless critical loads of Gerber beam. 


buckling load of dimensionless nonlocal Euler beam, Ppyr: critical buckling load of 
dimensionless nonlocal Timoshenko beam.) 
Dimensionless parameters are selected as follows: 


y b h 


g 


where b and h are the dimensions of the beam cross-section. 


5. Concluding Remarks 


In this work, based on Eringen’s nonlocal elasticity theory, the buckling of 
Timoshenko beam with the method of initial values is studied. The nonlocal 
Timoshenko theory considers both the scale effect and the effect of transverse shear 
deformation. For this reason, it has a significant role in the investigation of nanoscale 
elements. A priority of the method is that the high degree of statical indeterminacy 
adds no further hardship to the solution of the problem.” It is interesting to note 
that the size of the matrix that gives the buckling determinant is 2 x 2. If the 
conventional technique applies to these problems the buckling loads are found by the 
determinant of 4x4 matrix (see Sec. 4.1 and 4.2). For the problems given in 
Secs. 4.3 and 4.4, the buckling determinant was obtained with the help of 3 x 3 
determinants. 8 x 8 determinants arise in the solution of the same problems with the 
conventional technique. Nonlocal theories do not neglect small-scale effects, for this 
reason, the critical loads calculated by local and non-local theories are different for 
large values of y/L (see Figs. 4-12). In the nonlocal Timoshenko beam, both the 
scale effect and the shear effect are considered. For this reason, the safest critical 
loads are those belonging to the nonlocal Timoshenko beam (see Figs. 4-12). 
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